A  D- 759  549 


THEORETICAL  CALCULATION  OF  ACOUSTIC- 
GRAVITY  WAVE  EXCITATION  BY  NUCLEAR 
EXPLOSIONS  IN  THE  ATMOSPHERE 


C.  A.  Newton,  et  al 


Teledyne  Geotech 


Prepared  for : 

Air  Force  Office  of  Scientific  Research 
5  March  1973 


DISTRIBUTED  BY- 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


ad  759549 


AL-72-4 


. CMitrihtlif  ti  mtt'i 

n^nl»%  t!  It*  minmmtut  w »rM 


D  D“C\ 

[(JE0J3S.GE 


f] 

M4V  «  1973  j  I 

fE©E9r;;  uJIJJ 

THEORETICAL  CALCULATION  OF 
ACOUSTIC-GRAVITY  WAVE  EXCITATION 
BY  NUCLEAR  EXPLOSIONS  IN  THE 

ATMOSPHERE 


BY 

C.A.  NEWTON.  E.A.  FlINN  and  D  O  HARKRIOER 

SPONSORED  BY 
ADVANCED  RESEARCH  PROJECTS  AGENCY 

MONITORED  BY 
AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 

ARPA  ORDER  NO.  1357 
5  MARCH  1973 


K«pfoduc«a  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 
U  5  v ACrm*"* 


"WTELEDYNE  GEOTECH 

ALEXANDRIA  LABORATORIES 


Approved  for 
distribution 


Public  release, 

unlimited. 


w_ _ (S*tu:“y  clftlllcalton  ol  title,  b otfy  «<  nhetrect  end  IndmMln#  |-(>(»*0on  trutl  he  entered  »trn  the  overall  report  It  cl*  filled) 

I  OHIIi  NATIHC  ACTIVITY  (Corporate  author)  '  '  •  I*..  REMO*  ,  SFCUHITV  CLASSIFICATION - 

Tel* dyne  Gcotech  UNCLASSIFIED 

‘  Alexandria,  Virginia  22313  —— — - 


Ufc.  CROUP 


L  3  REPORT  TITLE 


THEORETICAL  CALCULATION  OF  ACOUSTIC -GRAVITY  WAVE  EXCITATION  BY  NUCLEAR  EXPLOSIONS 
IN  THE  ATMOSPHERE 


|  4  DESCRIPTIVE  NOTES  (Typ*  ot  rapirl  mnd  lnctu»l*m  dm  fa) 

|  Scientific  Interim 

fr  authoriii  (f'ltet  name,  middle  Initial,  latt  name) 

I  C.  A.  Newton 

E.  A.  Flinn 
D.  G.  Harkrider 


6  REPORT  DATE 

5  Mar  1973 

»  CONTRACT  OR  CHANT  NO 

F44620-69-C-0082 

6.  PROJECT  NO 

AO  1357 

6  270  ID 


10.  DISTRIBUTION  STATEMENT 


7«.  TOT 4L  NO.  OF  PACES  lb.  NO.  OF  REFS 

I  4  _ 

•a.  ORIGINATOR’*  REPORT  NUMBER,*) 

AL- 71- L/ 


t6.  OTHER  REPORT  NOH)  (Any  othmr  numbmrm  thmt  mmy  bm  ma*l0nmd 
thl§  report) 

iitm-XR-'  7  s  «  0  7  1  4 


Approved  for  Public  Release;  distribution  unlimited. 


Ill  SUPPLEMENTARY  NOTES 


Isa  SPONSORING  MILITASV  ACTIVITY 


ArusKfNt'r) 

|  TECH,  OTHER  1400  Wilson  Boulevard 

j  Arlington,  Virginia  22209 

[IS-  attract  - 

t  Jrom  theoretical  calculations  of  acoustic-gravity  waves  in  the  period  range  40  to 
1,000  seconds  from  various  combinations  of  yield  and  burst  height,  we  conclude  that 
the  power  varies  as  the  square  of  the  yield  for  low-altitude  explosions  and  as  the 
two-thirds  power  for  altitudes  above  100  km.  For  a  given  yield,  the  power  increases 
with  burst  height  to  a  maximum  which  is  found  at  lower  altitudes  for  higher  yields. 
The  affect  of  variable  atmospheric  structure  along  the  propagation  paths  render  it 
possible  to  make  only  a  rough  estimate  of  yield  and  no  meaningful  estimate  of  burst 
height.  The  study  also  shows  that  all  of  the  first  ten  modes  are  always  relatively 
well  excited  so  that  the  relative  modal  structure  of  a  signal  is  not  a  good 
diagnostic  for  burst  height. 

\ 


DD 


,1473 


UNCLASSIFIED 


Security  CUnRificalion 


Acoustic-Gravity  Waves 
Nuclear  Explosions 


THEORETICAL  CALCULATION  OF  ACOUSTIC-GRAVITY  WAVE 
EXCITATION  BY  NUCLEAR  EXPLOSIONS  IN  THE  ATMOSPHERE 

ALEXANDRIA  LABORATORIES  REPORT  NO.  AL-72-4 


Effective  Date  of  Contract: 
Contract  Expiration  Date: 
Amount  of  Contract  Dollars: 
Program  Code: 

Contract  Number: 

ARPA  Order  No.: 

Principal  Investigator: 
Period  Covered: 


1  March  1969 
31  March  1972 
$  590,572 
1F10 

F-44620-69-C-0032 

1357 

Carl  A.  Newton 

1  March  1969 
through 
31  March  1972 


Approved  for 
distribution 


Pub  He  release, 
unlimited. 


1C 


ABSTRACT 


Prom  theoretical  calculations  of  acoustic-gravity 
waves  in  the  period  range  40  to  1,000  seconds  from 
various  combinations  of  yield  and  burst  height,  we 
conclude  that  the  power  varies  as  the  square  of  the 
yield  for  low-altitude  explosions  and  as  the  two-thirds 
power  for  altitudes  above  100  km.  For  a  given  yield,  the 
power  increases  with  burst  height  to  a  maximum  which  is 
found  at  lower  altitudes  for  higher  yields.  The  affect 
of  variable  atmospheric  structure  along  the  propagation 
paths  render  it  possible  to  make  only  a  rough  estimate 
of  yield  and  no  meaningful  estimate  of  burst  height .  The 
study  also  shows  that  all  of  the  first  ten  modes  are 
always  relatively  well  excited  so  that  the  relative 
modal  structure  of  a  signal  is  not  a  good  diagnostic 
for  burst  height. 
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INTRODUCTION 


The  continuing  tests  of  nuclear  weapons  in  the 
atmosphere  have  prompted  interest  in  using  infrasonic 
data  recorded  at  distant  stations  to  estimate 
yields  and  burst  heights  of  the  sources.  We  have  calcu¬ 
lated  theoretical  barograms  for  a  variety  of  yields  and 
burst  heights,  synthesizing  the  barograms  from  the 
acoustic-gravity  wave  modes  which  propagate  in  the 
period  range  of  the  operating  microbaro^raph  arrays,  i.e., 
periods  of  4  to  14  minutes.  We  show  examples  of  barograms 
recorded  at  teleseismic  distances  for  many  combinations 
of  heights  and  yields,  and  draw  tentative  conclusions 
about  diagnostics  which  might  be  deduced  from  such 
recordings . 
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METHOD 


The  digital  computer  programs  used  were  those 
written  by  Harkrider  for  calculation  of  the  modal 
excitation  of  acoustic-gravity  waves  (Harkrider,  1964; 
Harkrider  and  Flinn,  1970).  The  original  programs  used 
a  mass - in j ect ion  source;  we  modified  these  programs  to 
use  instead  Pierce's  energy  injection  source  (Pierce, 
1968)  which  is  physically  more  realistic. 

Using  these  programs  it  was  straightforward  to 
calculate  the  spectrum  and  hence  the  time-domain  waveform 
(via  Aki's  synthesis  technique)  for  arbitrary  yields 
and  burst  heights,  for  any  of  the  various  acoustic  or 
gravity  modes. 

The  atmospheric  model  used  was  the  standard  ARDC 
model  (Wares  et  al.,  1960),  bounded  below  by  a  rigid 
surface  and  above  by  a  free  surface.  The  calculated 
phase  velocity  dispersion  curves  for  this  model  are 
shown  in  Figure  1. 

For  standard  instrumentation  in  use,  and  at  an 
epicentral  distance  of  10,000  km,  only  the  fundamental 
gravity  mode  GRq  and  the  lowest  acoustic  modes 
through  S3  are  involved  in  synthesis  of  waveforms.  We 
calculated  theoretical  barograms  cor  each  of  these 
gravity  and  acoustic  modes,  as  well  as  their  relative 
excitation  for  each  yield  and  source  height.  The 
composite  barograms  in  Figures  2  through  10  were 
constructed  from  the  scaled  modal  barograms.  The  notches 
in  the  spectra  are  due  to  the  wave-guide  propagation, 
but  we  believe  that  the  ripple  is  probably  an  artifact 
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of  the  synthesis  algorithm  we  employed. 

The  maxima  in  the  long-period  spectra  (1  minute) 
are  shown  in  Figures  11  through  13.  The  general  conclu¬ 
sions  we  draw  from  these  figures  are  that:  (1)  the  long- 
period  power  varies  as  the  square  of  the  yield  for  low- 
altitude  explosions,  and  as  the  two- thirds  power  for 
high  altitudes  (h>100  km);  (2)  for  a  given  yield,  the 
long-period  power  increases  with  burst  height  to  a 
maximum,  which  occurs  at  lower  altitudes  for  higher 
yields . 

Because  of  winds  and  nonuniform  atmospheric  struc¬ 
ture  along  the  signal  propagation  paths  (effects  which 
we  have  not  taken  into  account  in  the  present  work) 
only  a  rough  estimate  of  yield  can  be  made  from  the 
observational  data,  and  no  meaningful  estimate  of  burst 
height  can  be  attempted. 

Acoustic  spectrograms  from  atmospheric  explosions 
contain  energy  at  periods  as  low  as  a  few  seconds.  In 
order  to  identify  the  corresponding  acoustic-gravity 
wave  modes,  we  calculated  the  medium  response  for  all 
the  acoustic  modes  from  Sq  to  S^g,  at  periods  down  to 
a  few  seconds.  We  chose  ten  modes  to  include  in  a  synthe¬ 
sis  of  the  broadband  waveforms.  The  ten  modes  were  GRq 
and  Sq  through  Sg)  over  the  frequency  intervals  shown 
in  Figure  14.  To  test  the  importance  of  the  higher 
acoustic  modes  for  identifying  source  height,  we  calcu¬ 
lated  theoretical  barograms  for  five  source  altitudes 
from  1  km  to  45  km.  The  theoretical  barograms  for  each  of 
ten  modes  for  a  source  yield  of  IOC  kT  at  an  altitude  of 
5  km  are  shown  in  Figure  15.  The  moving-window  power 
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spectra  for  a  composite  barogram  consisting  of  5  modes  is 
shown  in  Figure  16;  Figure  17  shows  a  ten-mode  composite. 

The  composite  waveforms  for  five-mode  and  ten-mode 
synthesis  are  shown  in  Figures  18  through  20  for  yields 
of  1  kT  and  100  kT  and  altitudes  of  1  km,  5  km,  and  15  km. 
It  is  clear  from  observational  data  that  attenuation  and 
scattering  reduces  the  higher  mode  amplitudes  much  more 
than  the  lower  mode  orders.  The  high-frequency  ripple 
seen  in  the  ten-mode  composites  is  an  artifact  introduced 
by  the  relatively  narrow  frequency  interval  in  the 
synthesis. 

The  peak  overpressures  for  each  of  the  ten  modes  are 
shown  in  Figure  21.  This  shows  tnat  all  ten  modes  are 
relatively  well  excited,  and  makes  it  clear  that  relative 
modal  overpressure  is  not  a  good  diagnostic  for  burst 
height . 
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Figure  3.  Synthesized  barograns  and  spectra  for  sources  at  5  km 
altitude  and  yields  from  1  kT  to  50  MT.  Each  record  extends  over 
the  group  velocities  from  333  to  290  m/s,  beginning  400  minutes 
after  an  event  8000  km  distant.  The  horizontal  axes  are  scaled 
m  seconds  of  period. 


Fieurc  4.  Synthesized  barograms  and  spectra  for  sources  at  10  km 
altitude  and  yields  from  1  kT  to  50  MT.  Each  record  extendr  over 
the  group  velocities  from  555  to  290  m/s,  beginning  *00  »ctcrs 
after  an  event  8000  km  distant.  The  horizontal  axes  arc  sealed 

in  seconds  of  period. 


Figure  5.  Synthesized  barograas  and  spectra  for  sources  at  20  km 
altitude  and  yields  fro*  1  kT  to  50  MT.  Each  record  extends  over 
the  group  velocities  fro*  333  to  290  */s,  beginning  400  ainutes 
after  an  event  8000  k*  distant.  The  horizontal  axes  are  scaled  in 
seconds  of  period. 


11 


Yiw 

M 


hi'/ 1  '  1  \ 


,/  1  Kt 


W/7  M 


10  Kt 


f \J.  .  .  .  D  u  . 


: AH  V, 


100  Kt 


- 1  ...5  V 


.  .  v  *-  • 


1“~ 


10  Mt 


60  Mt 


Figure  6.  Synthesized  barogrw*  *^  MT.  '^h  ^(Tminutes 

ifes-2?Sffii&  -■ — 

in  seconds  of  period. 


*rj?r 


> 


100  Kt 


Figure  7.  Synthesized  barograms  and  spectra  for  sources  at  40  km 
altitude  and  yields  from  100  kT  to  50  MT.  Each  record  extends  over 
the  group  velocities  from  333  to  290  m/s,  beginning  400  minutes 
after  an  event  8000  km  distant.  The  horizontal  axes  are  scaled  in 
seconds  of  period. 
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Figure  8.  Synthesized  barograms  and  spectra  for  sources  at  50  km 
altitude  for  yields  from  100  kT  to  50  MT.  The  horizontal  axes  are 
scaled  in  seconds  of  period.  Each  record  extends  over  the  group 
velocities  from  333  to  290  m/s,  beginning  400  minutes  after  an 
event  8000  km  distant. 


Figure  9.  Synthesized  barograms  and  spectra  for  sources  at  75  km 
altitude  and  yields  from  1  kT  to  50  MT.  Each  record  extends  over 
the  group  velocities  from  333  to  290  ra/s,  beginning  400  minutes 
after  an  event  8000  km  distant.  The  horizontal  axes  are  scaled  in 
seconds  of  period. 
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Figure  10.  Synthesized  barograms  and  spectra  for  sources  at  100  km 
altitude  for  yields  from  1  kT  to  50  MT.  Each  record  extends  over 
the  group  velocities  from  333  to  290  m/s,  beginning  400  minutes 
after  an  event  8000  km  distant.  The  horizontal  axes  are  scaled  in 
seconds  of  period. 


111 

iWlr 

m 

mum. 

UKUf 

jiyife 

mm, 

iiMin 

FREQUENCY  (cpi) 


°f  acoys  tic -gravity  inodes  for  a  unit  source 
amplitude.  The  hatchured  regions  show  the  frequencies  where  modal 

1Sinr!ater  than  10-2  *  and  the  bars  indicate  excitation 
greater  than  10  °9 


trograa  for  the  synthetic  barograa 
attitude,  with  five  aodes. 


■odes 


synthesis 


v)  o 
E  a 

O 

«  a 

1-  3 

* 

4-) 

M 

O  V 

1- 

l-i  X 

o 

n  h 
J3 

i 

o 

•••til  akin 


modes 


ALTITUDE  (Km) 


Figure  23.  Variation  of  peak  amplitude  with  altitude  for  the 
individual  modes  used  in  the  barogram  synthesis. 


